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Abstract 

In this paper we revisit and extend a previous analysis where the possible relevance of 
quantum gravity effects in a cosmological setup was studied. The object of interest are non- 
| local (logarithmic) terms generated in the effective action of gravity due to the exchange 

f^) . in loops of masslcss modes (such as photons or the gravitons themselves). We correct 

one mistake existing in the previous work and discuss the issue in a more general setting 
in different cosmological scenarios. We obtain the onedoop quantum-corrected evolution 
equations for the cosmological scale factor up to a given order in a derivative expansion in 
two particular cases: a matter dominated universe with vanishing cosmological constant, 
CO ' and in a de Sitter universe. We show that the quantum corrections, albeit tiny, may have 

a secular effect that eventually modifies the expansion rate. For a de Sitter universe they 
(j ' tend to slow down the rate of the expansion, while the effect seems to have the opposite 

O" 1 sign in a matter dominated universe. To partly understand these effects we provide a 

5^ ' complementary newtonian analysis. 

wo: 

1 Introduction 

> 

It has been said pQ that the effective action of quantum gravity is the most effective of all 
effective actions, meaning that an expansion in powers of p 2 /16nMp would give in normal 
conditions such a tiny contribution to any scattering amplitude that the 0{p 2 ) (the usual 
Einstein-Hilbert) term is good enough for all practical purposes (and even for many non- 
practical ones). This is unlike pion physics where the presence of higher order operators leads 
lo measurable effects already at moderate energies. Thus, the fact that the 0(p 4 ) terms are 
ultraviolet divergent does not really bother us in practical calculation^, although of course 
^ . the issue is very relevant from a fundamental point of view. 

It is easy to see why quantum corrections are so small. In fact, as already implied above, 
the expansion is on powers of p 2 /l6ir Mp (actually V 2 /lQirMp, 1Z/16ttM p ) and therefore very 
small for physical values of the energy or curvature. Non local pieces in the effective action 
(~ In V 2 ), due to the presence of strictly massless modes, somewhat increase the relevance of 
higher order terms, but locally they are still negligible. 

There are two reasons why such apparently hopelessly small corrections might nevertheless 
be relevant in a cosmological setting. One reason is that curvature was much larger at early 
stages of the universe. For instance, in a de Sitter universe 1Z ~ H 2 . In a inflationary scenario 
H 2 = 87rGVo/3, Vq being the scale of inflation that is limited by CMB measurements to be 
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1 In pion physics, and in quantum gravity too, one can make sense of non-renormalizable interactions at the 
expense of introducing more subtraction constants, one per independent operator in the derivative expansion. 
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V < 10 16 GeV. The Hubble constant H could have been as large as H ~ 10 13 GeV (the 
present value is 10 2 GeV). However, even in the most favourable case, the correction is still 
nominally of (D(10~ 12 ) or less and it seems too small to be seen. Or maybe not? 

Indeed a second reason to study this problem is that this nominal suppression overlooks 
the fact that the logarithmic non local term corresponds to an interaction between geometries 
that is long-range in time, an effect that does not have an easy classical interpretation. When 
integrated over time, this may bring about a large enhancement of the contribution of higher 
order contributions to the point where new interesting affects may appear. 

It is quite important to realize that there is no ambiguity in the overall coefficient of 
the logarithmic non-local term as it depends only on the structure of the Einstein-Hilbert 
Lagrangian and the number of massless modes (or modes whose mass is much smaller than 
the inverse of the horizon radius). Thus the effects are model independent and, if observable, 
can be unambigously predicted, at least inasmuch as one can make accurate predictions in an 
effective theory. 

The possibility of observable effects of the non-local terms in the effective action in such a 
situation was recognized in [2]. Although the results in [2] appear to be correct in substance, an 
error slipped in the calculation, unfortunately. This error is corrected here. More importantly 
the analysis is extended to different cosmological models. It is seen that a secular effect from 
the non-local terms can be predicted unambigously (within the approximations inherent to an 
effective action, that is up to terms with higher derivatives) and it is seen to lead to potentially 
visible effects. 

The relevance of quantum gravity corrections on inflation was also pointed out in [3]. 
Unfortunately, it is difficult to draw a parallelism between the two approaches. For one thing, 
we are finding here a one- loop effect, while the one discussed in [3] is a two-loop one due to 
particle creation [1] and thus clearly subleading. 

The importance of non-local term^l in the effective action of gravity cannot be overem- 
phasized. This has been recently reviewed in [5], although the non-localities discussed in that 
paper do not actually correspond to the one present here. 

In section [2] we rederive the loop-corrected evolution equation for the cosmological scale 
factor in the presence of non-local logarithmic terms. In section [3] we apply these techniques 
to a matter dominated universe, governed by a power-law expansion. In section [3] we rederive 
the quantum corrections to the cosmological evolution equation in a de Sitter background. In 
section we present the numerical analysis of the solutions and comment on their physical 
relevance. In order to better understand the relevance of the logarithmic terms we have 
provided a classical analogy that is presented in section El 

The metric convention we use in this paper for Minkowski space is 





The Einstein equations are 




8irGT, 



A 9, 



(2) 



2 By non-local we mean terms non-analytic in V 2 , such as the In V 2 pieces that appear in one-loop effective 
gravity. We do not consider the so-called f(7V) gravity. 
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where g pv is the metric tensor, A is the cosmological constant, R = g^ u R pu and 

p a pa _ o -pa , -pa p/3 _ pa p/3 /o\ 

T l/3 = \a 1P ( d f39pa + d a g p/3 - d p g a p) . (4) 
The previous equations are derived from the action 

S = j dx "J-g{K - 2A) + S ma tter- (5) 

Quantum corrections to the Einstein-Hilbert action were originally computed by 't Hooft 
and Veltman in [B] in the case of vanishing cosmological constant, and by Chistensen and 
Duff for a de Sitter background [7]. Other related references that we have found particularly 
useful in the present context are [8] and [9]. 

The key ingredient we shall need is the divergent part of the one-loop effective action. 
Using dimensional regularization and setting d = 4 + 2e we get [9] 

rf f v f = J dxV^giciK^W + c 2 A 2 + c 3 TZA + c 4 ft 2 ]. (6) 

The constants Cj are actually gauge dependent and only a combination of them is gauge 
invariant. Using the equations of motion (in absence of matter) TZ^ V = g^ v A, the previous 
equation reduces to the (gauge-invariant) on-shell expression 



rS^/^fA*. (7) 



If we set A = in ©, we get the well-known 't Hooft and Veltman divergence, that in the 
so-called minimal gauge is 

Yi% = V I dx^=g[—K uv K t * v + — K 2 }. (8) 

e// 16vr 2 e J v L 20 M 120 J K J 

If the equations of motion are used in the absence of matter this divergence is absent, as is 
well known. 

Once the divergence is determined it is straightforward to obtain the logarithmic pieces 
since they always appear in the combination 

1 V 2 

- + ln-V. (9) 

The dimensionful constant [i is the substraction scale that is required for dimensional consis- 
tency. 



2 The equations of motion in the presence of non-local terms 

In this section we shall derive quantum corrected equations of motion for the cosmological scale 
factor including the non-local logarithmic terms that appear in the one-loop effective action. 
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We shall consider here a simplified effective action that includes only terms containing the 
scalar curvature we split the action into three parts and redefine the constants for convenience 

S = k 2 f j dx^/^gTZ + a J dx^/^gTZln(V 2 / fi 2 )TZ + (3 J dx^/^gTZ 2 ^ 

= k 2 (Sx + aS 2 + pS 3 ) , (10) 



where n 2 = Mp/IQ-K = 1/lQirG. \i is the subtraction scale whose contribution is by itself local, 
but gives the right dimensions to the non-local term. The coupling /3 is /i dependent in such 
a way that the total action S is /i- independent. While the value of (3 is actually dependent 
on the UV structure of the theory (it contains information on all the modes -massive or not- 
that have been integrated out), the value of a is unambiguous as it depends only on the IR 
structure of gravity, that is entirely described by the Einstein-Hilbert Lagrangian and the 
massless modes. 

In conformal time, dt = adr, we have 



9p,v = a (t)^, K = Q^tt-t, \f ZI 9 = a (r). (11) 

In order to obtain the modified equations of motion for the cosmological scale factor, we 
first perform the variation of the local action, getting the following results 

6S, 



5a(r 



\ = 12a" (12) 



5S 3 / (a") 2 a' a!" ia'fa" a< 4 > 



72 - - 4-^ + 6^j- + -j- (13) 



Sa(r) I a 3 a 3 a 4 a 

The d'Alembertian in conformal space is related to the Minkowski space operator by 

V 2 = cT 3 Da + -ft (14) 
6 

Neglecting the 1Z term in the previous equation and commuting the scale factor a with the flat 
d'Alembertian (terms with higher derivatives are neglected in the effective action philosophy), 
we can write 

v2 ^r n _ <i5> 

Where ao = a(0). With this rescaling (absorbable in (3), at r = the d'Alembertian in 
conformal space matches with the Minkowskian one. We will set ao = 1 for simplicity. 

We can now separate S 2 in turn into a local and a genuinely non-local piece 
S 2 = J dx^j {-2n\n(a)lZ + lZ\n{U/iJL 2 )U) 

^Si + Si 1 . (16) 

The variation of 5^ gives 

- - 72 ( ^f- [12 ln(a) - 10] + ^ [-8 ln(a) + 4] + (17) 



6*a(r 



(a") 2 

[_ 6 ln(a) + 2] + -=-2 ln(a) } (18) 
a z 
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In order to determine the variation of S^ 1 we need to compute the Green function 



{x\ \n\3\y). 



(19) 



We follow the method of [2] that we shall not repeat here. We mention here that the normal- 
ization of the delta function used in [2] was non-covarianlH. Using the proper normalization 
and correcting for this mistake we find 



Si 1 = 36 I dr 



a"{r) [* 1 a»(r>) 



a i T ) Jo t-t' a(r') 



(20) 



The variation of Si 1 is 



5S. 



ii 



5a{r) 



36 



2a~ 3 (r) (a / (r)) 2 -2a- 2 (r)a"(r) 



dr 



, 1 a'V) 
r — r' a(r') 



-2a-\T)a!(T) 



1, n 9 ( r , / 1 a "( T ')\ -i/ n 92 ( r , / 1 a "( r 

dr 



Ot 



dr' ' , V +a" i (r) 
t — r a(T n 



8t 2 



t — t' a(r') 



(21) 



Using repeatedly integration by parts it is possible to further simplify the expression elimi- 
nating the derivatives acting on the integrals 



+ / dr 



d_ 

dr 

T ^ 1 



d.T 



, 1 a"(r') 



-i„// 



-+ 



r — r 



r — r' a(r') 
7 [-a-V)«V)aV) + a-VVV)] 



(22) 



9r 2 \7 r — r' a(r') / r 2 

+ [dr*—, [2a-\r'){a'{r')fa"(r')-a-\r'){a'\r')f 
Jo r — r 

-2a- 2 (r')a'(r')a'"(r>) + a -\r')a^ (r')\ (23) 

The terms that potentially diverge at r = arise from the fact that we patch together 
Minkowski and FRW spaces at r = 0. If it is done smoothly enough, the derivatives of the 
scale factor should vanish at that point. 

The previous expressions are formal in the sense that the short distance singularities 
have not been properly regularized. To regulate the region where r — ► r' we use dimensional 
regularization as explained in detail in [2] . The regulator is in fact introduced when the Green 
function (|19|) is computed in 4 + 2e dimensions. In practice, this amounts to the replacement 

1 u~ 2e 

(24) 



( T _ T /)l+26> 



where fi is the subtraction scale previously introduced. The physical result is obtained for 
e — > 0. The logarithm of the effective action (jlOf) is then reproduced and the 1/e divergence 
(proportional to 1Z 2 ) is cancelled by the (divergent) counterm included in 0. Thus it is easy 



3 We thank G.Perez for bringing this to our attention 
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to compute, even numerically, the genuinely non-local piece, as the would-be singular term in 
the integral is clearly identified. 

At this point there are several ways to proceed. One might of course attempt to find 
solutions of the equation of motion for a(r) obtained by adding the variations for Si, /3S*3 and 
0(5*2 + that we have just computed. This way of proceeding is not really justified if the 
0(p 4 ) terms are understood correction. 

In the spirit of effective Lagrangians it is better to proceed otherwise. We obtain first the 
lowest order equation of motion from Si and plug it in ^(Sg + S| ) + The quantum 

corrections act then as an external driving force superimposed to Einstein equations. This 
procedure of course gives trivially a net zero additional contribution in the present toy model 
as neither matter nor a cosmological constant have been considered. In the next sections we 
shall introduce T^ IU and A to find more interesting effects. 



3 Quantum gravity effects in a matter dominated universe 

Let us consider a pressureless distribution of matter characterized by the energy-momentum 
tensor T^ u = diag(p, 0, 0, 0). We can use Einstein equations to write (in conformal time) 

TZ 00 = -4vrGa 2 (r)p, fty = 4vrGa 2 (r)p. (25) 

The addition of an energy-momentum tensor leads in principle to the appearence of 1ZT, T 2 
and T^yT^ terms in the effective action. However, the details of the calculation are still 
controversial to some extent. Here we shall adopt the view that a classical energy- momentum 
tensor (such as one representing presureless dust) interacts only via classical gravitons. 

We then substitute these expressions in the non-local 0{p A ) term effective action 

S = — — / dxy/^TZ (26) 



16ttG 
1 



16vr 



dx^—g 



7 V 2 IV 2 

— W v In -fH",, + n In -^rll 

20 v fi 2 120 fi 2 

A 



(27) 



+ local terms of 0(p ) + S ma tter- (28) 

The local 0(p 4 ) terms (proportional to V? and lV' 1 /R, h ' ^ can be used to absorb the divergences 
of the one-loop action. If one day one would be able to make a precise measurement of some 
quantum gravity effect, one could determine these coefficients and all other predictions would 
be calculable unambiguously with an C(p 4 ) precision. For the time being, we only know 
unambigously the logarithmic non-local terms. But these are the ones that are like to give 
enhanced contributions as we discussed in the introduction. 

After use of the lowest order, 0(p 2 ) equations of motion, the non local terms simplify 
considerably 

S = k 2 (J dx^gU. + a J dx^gpln^-p^ + S mat ter (29) 
= k 2 ( dx^gTZ + aS 2 ) + S mat ter- (30) 



We omit the purely local part (analogous to S3 in the previous section) as is not calculable 
from the low energy information only and logs dominate anyway. 
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The value of a can be determined readily from the results of 't Hooft and Veltmann [6] 
(which includes only the contribution of virtual gravitons) and the lowest order equations of 
motion in the presence of matter 

a = -WttG 3 x — . (31) 
30 V ; 

The correction from massless photons or Yang-Mills fields does not seem to change the sign 
of a |10j . We have not considered other possibilities. In fact, the precise value of a is not so 
important (but the sign and its rough order of magnitude is). 

S2 modifies the equations of motion by adding a new term that is simply 



£py = -2p(0 2 a 3 (r)[41n(a(r)) + l] + 2p(r)a(T) J p(r')a 2 (/) » e , 



multiplied by a. Notice that quantum effects introduce long range interactions in time between 
global matter densities at different times. 

The value of p{r) is known from the lowest order equation of motion. In conformal time 

p(r) ~ a" 3 (r), a(r)~r 2 . (33) 

We shall discuss the physical relevance of these corrections after discussing in detail the 
solution in the de Sitter case. 



4 Quantum gravity effects in a de Sitter universe 

We shall proceed in a way similar to the previous section, but we shall now omit the energy- 
momentum tensor in the lowest order equations of motion. In fact it is not correct to simply 
assume T^ u ~ g^A in the equations of motion and just use the previous formulae (that would 
roughly be equivalent to exchanging p by A in the previous section, up to constants). The 
cosmological constant is the most relevant operator in gravity and it must be introduced from 
the outset. 

The relevant one-loop corrected effective action is 

\ f If 29 

S = / dxy/=g(K - 2A) + =■ / dx^/^g— A In -5- A + local terms of C(p 4 ) (34) 

167rG J I67H J 5 p z 

= k 2 ( / dx^f=g(K- 2A) + aS 2 J . (35) 



Now a is very different from the previous case 

G 29 

a = - x — (36) 

7T 5 

The dimensions of a are of course different as the dimensionality of p and A is not the same. 
Most importantly, it has the opposite sign. Notice that all our expressions are written in 
such a way that it is possible to consider a time-dependent cosmological constant (or matter 
density). 
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We split 52 in two parts 



S : 2 =-2 j dx^^k 2 \n{a) (37) 
S 1 ^ = J dxy^gAln(D/n 2 )A, (38) 



and obtain the corresponding variations following the method outlined in section [2] 



2 - -2A 2 a 3 (r)[41n(a(r)) + l] (39) 



<5a(r) 



6S " - 2A 2 a(r) f dr'a\r') r _ !^— . (40) 



Sa(r) 7 \ T ~ t 

The equation of motion will be 

12a" (r) - 8Aa 3 (r) + d-^ 2 - = 0, (41) 

oa{T) 

which at lowest order is just 

12a" (r) - 24H 2 ci 3 {t) = 0, (42) 
where H 2 = A/3. The solution of flM} is 

a:(r) _ ^ (43) 

The final step to solve iteratively ([IT]) is to plug the 0-th order solution aj(r) into the variation 
of 52 and recalculate the solution for a(r). 

It is clear that, apart from the sign difference, the quantum effects are formally very 
similar for a matter dominated and for a de Sitter universe. However, the fact that the signs 
are opposite means that their back reaction is completely opposite. If quantum corrections 
enhance expansion in one case, they will slow it down in the other. Futhermore, the size of 
the correctios is very different: the corrections in a matter dominated universe are down by 
a factor H 2 /Mp with respect to the ones in a de Sitter space-time with a large cosmological 
constant, already expected to be small. Let us now investigate the numerical relevance of the 
latter ones. 



5 Solving the evolution equation 

As we just discussed, we proceed by solving the varied gravitational action by a perturbative 
approximation, i.e., we consider the non-standard terms as a correction to the standard infla- 
tionary solution. This perturbative procedure is of course only valid as long as the correction 
is small compared to the unperturbed solutions. 

Before doing that we find it convenient to change time coordinates by introducing a variable 
s defined through aj = e s . Then 

= -2A 2 e 3s [4s + 1] (44) 

a/ 



5a{r) 
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5Sl_ 
8a(r) 



2AV/(s) 



a>i 



and the equation of motion reads 

e 2s a"(s) + e 2s a'(s) - 2o 3 (s) = ^aH 2 (-e 3s (l + 4s) + e s I(s)) 

where I is defined in conformal time as 

aj(r') 



I(r) = ^ [ dr'- 
Jo 



T _ r /)l+2e 



1 

2e 



(r^)- 2e 2 ^i(l,2;l-2e;FT) 



(45) 



(46) 



(47) 



with 2-^1 being a hyper geometric function. Let us expand this expression around e = to 
the first order, disregarding higher orders since we eventually take the limit e — ► 0. Using 
2 F 1 (l,n; 1;Ht) = (1 - Hr)' n = oJ(r), we get 



/(r) =ln(r/x) 2^1(1,2;!;^) 



2 de 



2 Fi(l,2;l-2e;ffr) 



e=0 



ln(r/i)af(r) 



2 <9e 



2 F 1 (l,2;l-2e;^r) 



e=0 



This can be computed and written in s time as 



= ln ( |(l- e -s)) e ^ + e s (l 



se s ), 



and the equation to solve is 



e 2s a"{s) + e 2 V(s) - 2a 3 (s) = -aR 2 



(5s + 2)e is + e 2s + e 3s ln(-^(l 



(48) 



(49) 



(50) 



Note that a appears only in the combination aH 2 . Since there are H large uncertainties in 
H in practice only the sign of a is relevant. In addition, there is some ambiguity associated 
to the choice of the renormalization scale that appears in the combination \n(fi/H) 0. The 
dependence on the subtraction scale is mild (logarithmic) but it is inherent to the effective 
action approach. To reverse the sign of the effect one has to go to absolutely unreasonable 
values of fi. 

Eq. (I50p can now be solved for different values of H and fx, where H < 10 13 GeV. The 
solution is shown in Fig. [TJ We can see that the curves present a very similar behaviour for 
the different values shown, though a higher value of H leads earlier to deviations from the 
usual inflationary expansion. Higher values of (i also have this effect, which is larger as H 
increases (note that what is relevant in the equation is the ratio n/H). In fact, we can see 
from ([50]) that if we considered values of fj,/H large enough (but not relevant physically), the 
logarithm term would become dominant and the deviation would be positive. 



6 Newtonian approximation 

It is well known that one can derive the Friedmann equations using only newtonian physics 
[11]. We will redo this exercise considering quantum corrections to the Newton potential. 
These have been computed by several authors [12]. In fact, the magnitude of this correction 

4 Recall that a change of /i is equivalent to a redefinition of the local C(p 4 ) counterterms 
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H=10 13 H=10" H=10 9 



0.95 



0.9 




Figure 1: The scale factor relative to the inflationary expansion for different values of the 
renormalization scale p and the Hubble constant H (all units are GeV). s gives roughly the 
number of e-folds. The precise definition is given in the text. 




Figure 2: In order to derive the Friedmann equation, we consider a test particle on the surface 
of a virtual sphere in our infinite, homogeneous and isotropic universe. 

and even the sign has been the subject of a long controversy. The correct value for the quantum 
correction to the gravitational force between two masses is given for instance in [13] . This we 
believe to be the relevant modification in the context of this newtonian approximation, but 
we shall redo the analysis for an arbitrary value of the correction for the sake of completeness. 

Quantum corrections can be summarized in a modification of Newton's 'constant' in the 
following way 

G «= G ( 1 + ^)> ( 51 ) 
where £ is a dimensionless constant. 

We consider an infinite, homogeneous and isotropic expanding universe. We will describe 
the equations of movement for a test particle on the surface of a comoving sphere, expanding 
with the universe. The mass inside the sphere is p(4ir/3)r 3 , where p is the density of the 
universe and r the radius of the sphere, (see Fig. The kinetic energy of the (unit mass) 
particle is (l/2)r 2 . In order to compute its potential energy U we shall consider first the force 
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between our test particle and another particle with mass mat a distance I from it 



F = -G(r)m^ = -Gm^ - £G 2 hm-^ = F x + F 2 (52) 

The first part F\ is of course elementary and once integrated it leads to a contribution to the 
total potential energy U = U\ + U2 that is easy to guess 

GM 4nGp 2 
U 1 = — = — r (53) 

For the second part, we have to integrate F2 over the whole sphere and also the space outside 
of it, since we can not make use of Gauss theorem. After doing this we get the following result 
for the contributions from matter outside and inside the sphere, respectively 

F 2out = f (t&h) p 



(54) 



Where r is the unit radial vector pointing to increasing values of r. If we sum these two 
contributions, the divergences in the boundary cancel and we have 



U2 = -(^h)^ = -^(^n) p - (55) 



3r 3 3 vs ' 3 
This is a manifestation of Birkhoff 's theorem. 

We can now write the total energy of our test particle, that is 

E = -r 2 — ^r 2 - % {tG 2 K) P - (56) 
2 3 3 v ^ 1 3 v ; 

Writing r(t) = a(t)x and dividing by a 2 /2 we get the modified Friedmann equation. 

m 2 _ « + £ ((G 2 h] _ « ( 5T) 

a(t) J 3 H 3 ^ ; 3a 2 (t)x 2 a?{t) V ; 

Where K = —2E/x 2 is the curvature, which we set to since we are considering a flat space. 
As we see, the local coordinate x appears explicitly and in order to put precise numbers we 
have to be very precise as we define our unit system. This is not very relevant for practical 
purposes due to the smallness of the correction and we shall embed this into £. 

After identifying the cosmological constant A = 8irGp, we have 

iW y a a «o« 1 y (58) 



a(t)J 3 V 3 a 2 (t) 

In order to compare with the time coordinate we have introduced before, we revert to the 
variable s, remembering that, to the lowest order, H 2 = A/3. 

(aWe* = « 4 (*)(l + ^^) (59) 
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The resulting evolution can be computed easily and the effects are seen to be very small; 
so small in fact that to make them visible in a reasonable time evolution (comparable to the 
times seen in Fig. [T]) we have taken £Gh ~ 10~ 20 , which is an unreasonable value for sure 
(see Fig. [3|). The effect is to increase the expansion rate. This simple exercise teaches us 
two things. One is that it is incorrect to simple take the cosmological constant as part of the 
'matter balance' of the universe and stick to newtonian physics. Incidentally, the calculation 
would be the same for a matter dominated universe and the effects even tinier. The second 
thing we learn is that classical physics with a potential modified by quantum corrections is 
unable to reproduce the enhancement that the non-local terms bring about due to long range 
time correlations. 

7 Conclusions 

In this paper we have analyzed the relevance of the non-local quantum corrections due to 
the virtual exchange of gravitons and other massless modes to the evolution of the cosmo- 
logical scale factor in FRW universes. We have considered two different setups: a matter 
dominated universe, characterized by a matter density p(t), and a de Sitter universe with a 
large cosmological constant. 

In the de Sitter universe, while the effects are locally absolutely tiny, even after allowing 
for the largest possible value of H, we have found that they lead to a noticeable secular effect 
that slows down the inflationay expansion after a long time. This is a pure one-loop effect 
that is not actually related to particle creation and its back-reaction on the universe expansion 
and which constitutes a two loop effect as emphasized by Tsamis and Woodard. In a matter 
dominated universe the effect is a lot smaller, and it may be of the opposite sign. It is quite 
interesting that quantum effects seem to enhance the expansion rate in this case. 

The physical effects thus depend crucially on the sign of the quantum corrections. We 
have also seen that this effect has no really classical analogy. 

It is quite important to emphasize once more that the results presented here are not 'just 
another model'. Quantum gravity non-local loop corrections are required by unitarity, even 
if the theory is non-renormalizable; they can be computed quite precisely in a derivative 
expansion; and they appear to be of some relevance in the present situation. Perhaps the 
most interesting result of our analysis is indeed the fact that these effects can be predicted 
unambiguously within the limits of an effective theory. 
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